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Abstract 

An extensive Quantum Monte Carlo calculation is performed for the two-leg 
Hubbard ladder model to clarify whether the singlet pairing correlation decays 
slowly, which is predicted from the weak-coupling theory but controversial 
from numerical studies. Our result suggests that the discreteness of energy 
levels in finite systems affects the correlation enormously, where the enhanced 
pairing correlation is indeed detected if we make the energy levels of the 
bonding and anti-bonding bands lie close to each other at the Fermi level to 
mimic the thermodynamic limit. 
PACS numbers: 74.20.Mn, 71.10.Fd 
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Over the past several years, strongly correlated electrons on ladders have received much 
attention both theoretically and experimentallyS. This has been kicked off by the theoretical 
studies suggesting a formation of the spin gap and possible occurrence of superconductivity 
in such system Jii. 

The 

weak-coupling theory with the bosonization and renormalization-group techniques^Hi has 
indeed shown that the Hubbard model on a two-leg ladder has a spin gap, and if the system 
is free from umklapp processes, singlet pairing correlation function decays as ~ 1/r" with 
a = 1/2 (r: real space distance) in the weak-coupling limit. 

Since SDW and 2kp CDW correlations have to decay exponentially in the presence of a 
spin gap in a two-leg ladder, the only phase competing with superconductivity will be ikp 
CDW, whose correlation should decay as Hence the pairing correlation dominates 

over all the others if a < 1. 

As for the opening of the spin gap, the density-matrix renormalization group (DMRG) 
studies in the strong-coupling regime also indicate its presence in both t-J and Hubbard 
ladder mo delsSii. If we further focus on the t-J ladder, detects a pairing cor- 

relation decaying slightly slower than ~ 1/r and a CDW correlation decaying faster than 
~ 1/r for an electron density of n = 0.8 with J/t = 10. 

However, the dominance of the pairing correlation in the Hubbard ladder model seems to 
be a subtle problem in numerical calculations. Namely, a DMRG study by Noack et al. for 
the doped Hubbard ladder with n = 0.875, U/t = 8, and t± = t (where t and t±_ are intra- 
and interchain hoppings, respectively) shows no enhancement of the pairing correlation over 
the U = resultS, while they do find an enhancement at t± = Asai performed a 

Quantum Monte Carlo (QMC) calculation for a 36-rung ladder with n = 0.833, U/t = 2 and 
in which no enhancement of the pairing correlation was found. On the other 
hand, Yamaji et al. have found an enhancement for the values of the parameters where 
the lowest anti-bonding band levels for U = approaches the highest occupied bonding 
band levels, although their results have not been conclusive due to small system sizes (< 6 



runi 



.gs)Ei. 



Thus, existing analytical and numerical results appear to be controversial. This is dis- 
turbing since the superconductivity in the Hubbard ladder, especially with t± ~ t, is of great 
interest as a model for cuprate ladder-like materials, for which an occurrence of supercon- 
ductivity has indeed been reported very recentlyEl. In the present work, we have performed 
an extensive QMC calculation for the Hubbard ladder with ~ t in order to clarify the 
origin of the discrepancies among existing results. We conclude that the discreteness of en- 
ergy levels in finite systems affects the pairing correlation enormously, where the enhanced 
pairing correlation is indeed detected if we tune the parameters so as to align the discrete 
energy levels of bonding and anti-bonding bands at the Fermi level in order to mimic the 
thermodynamic limit. 

The Hamiltonian of the two-leg Hubbard ladder is given in standard notations as 



where a(= 1,2) specifies the chains. 

In the weak-coupling theory, the amplitude of the pair hopping process between the 
bonding and anti-bonding bands in momentum space flows into the strong-coupling regime 
upon renormalization, resulting in a formation of gaps in both of the two spin modes and a 
gap in one of the charge modes when the umklapp processes are irrelevant. This leaves one 
charge mode massless, where the mode is characterized by a critical exponent Kp, which 
should be close to unity in the weak-coupling regime. Then the correlation function of an 
interchain singlet pairing, Oj = (cij|C2i| — Cij|C2i|)/v^, decays like l/r^/('^^p\ 

Here, we have applied the projector Monte Carlo method^ to look into the ground 
state correlation function P(r) = {0\^^0i) of this pairing. We assume periodic boundary 
conditions along the chain direction, Cn+i = Ci, where is the number of rungs. 

The details of the QMC calculation are the following. We took the non-interacting Fermi 
sea as the trial state. The projection imaginary time r was taken to be ~ 60/t. We need 



n = -ty{c\ 




(1) 
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such a large r to ensure the convergence of especially the long-range part of the pairing 
correlation. This sharply contrasts with the situation for single chains, where r ~ 20/i 
suffices for the same sample length considered here. The large value of r, along with a 
large on-site repulsion U, makes the negative-sign problem serious, so that the calculation is 
feasible for U/t < 2. In the Trotter decomposition, the imaginary time increment [r/ (number 
of Trotter slices)] is taken to be < 0.1. We have concentrated on band fillings for which the 
closed-shell condition (no degeneracy in the non- interacting Fermi sea) is met. We set t = 1 
hereafter. 

We first show in Fig.|l]the result for P(r) for t± = 0.98 and t± = 1.03 with U = 1 and the 
band filling n = 0.867 = 52 electrons/ (30 rungs x 2 sites). The U = result (dashed line) 
for these two values of t± are identical because the Fermi sea remains unchanged. However, 
if we turn on U, the 5% change in t± = 0.98 — > 1.03 is enough to cause a dramatic change in 
the pairing correlation: the t± = 0.98 result has a large enhancement over the U = result 
at large distances, while the enhancement is not seen for t± = 1.03. 

In fact we have deliberately chosen these values to control the alignment of the discrete 
energy levels at f/ = in finite, two-band systems. Namely, when t± = 0.98, the one- 
electron energy levels of the bonding and anti-bonding bands for U = lie close to each 
other around the Fermi level with the level offset {Ae in the inset of Fig.|l|) being as small 
as 0.004, while they are staggered for t± = 1.03 with the level offset of 0.1. On the other 
hand, the size of the spin gap is known to be around 0.05t for U = 80, and is expected to be 
of the same order of magnitude or smaller for smaller values of U. The present result then 
suggests that if the level offset Ae is too large compared to the spin gap, the enhancement 
of the pairing correlation cannot be seen. By contrast, for a small enough Ae, by which an 
infinite system is mimicked, the enhancement is indeed detected as expected from the weak 
coupling theory, in which the spin gap is assumed to be infinitely large at the fixed point of 
the renormalization fiow. 

Our result is reminiscent of those obtained by Yamaji et al^, who found an enhancement 
of the pairing correlation in a restricted parameter regime where the lowest anti-bonding 
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levels approaches the highest occupied bonding levels. They conclude that superconductivity 
occurs when the anti-bonding band 'slightly touches' the Fermi level. However, our result 
in Fig.^ is obtained for the band filling for which no less than seven out of 30 anti-bonding 
levels are occupied at f/ = 0. Hence the enhancement of the pairing correlation is not 
restricted to the situation where the anti-bonding band edge touches the Fermi level. 

Now, let us more closely look into the form of P(r) for t± = 0.98. It is difficult to 
determine the decay exponent of P{r), but here we attempt to fit the data by assuming 
a trial function expected from the weak-coupling theory. Namely, we have fitted the data 
with the form P(r) = ^ J2d=±{'^^d + (2 — c)rj^ — [cos{2kprd) + cos{2kprd)] rj^} with 
the least-square fit (by taking logarithm of the data) c = 0.11. Because of the periodic 
boundary condition, we have to consider contributions from both ways around, so there are 
two distances between the 0-th and the r-th rung, i.e, r_|_ = r and r_ = N — r. The period of 
the cosine terms is assumed to be the non-interacting Fermi wave numbers of the bonding 
and the anti-bonding bands in analogy with the single-chain case. The overall decay should 
be 1/r^ as in the pure ID case. We have assumed the form c/r^/^ as the dominant part 
of the correlation at large distances because this is what is expected in the weak-coupling 
theory. A finite [/ ~ 1 may give some correction, but the result (solid line in Fig.|I|) fits 
to the numerical result surprisingly accurately. If we least-square fit the exponent itself as 
1/r", we have 0.2 < a < 0.7 with a similar accuracy. Thus a finite U may change a, but 
a > 1 may be excluded. To fit the short-range part of the data, we require non-oscillating 
(2 — c)/r^ term, which is not present in the weak-coupling theory. We believe that this is 
because the weak-coupling theory only concerns with the asymptotic form of the correlation 
functions. 

In Fig.^ we show a result for a larger system size (42 rungs) for a slightly different 
electron density, n = 0.905 with 76 electrons and tj_ = 0.99. We have again an excellent fit 
with c = 0.07 this time. 

In Fig.^, we display the result for a larger U = 2. We again have a long-ranged P{r) 
at large distances, although P{r) is slightly reduced from the result for U = 1. This is 
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consistent with the weak-couphng theory, in which Kp is a decreasing function of U so that 
once the spin gap opens for t/ > 0, the paring correlation decays faster for larger values of 
U. 

To explore the effect of umklapp processes, we now turn to the filling dependence for a 
fixed interaction U = 2. We have tuned the value of t± to ensure that the level offset (Ae) 
at the Fermi level is as small as 0(0. Olt) for f/ = 0. In this way, we can single out the effect 
of umklapp processes from those due to large values of Ae. If we first look at the half-filling 
(Fig-Ka)), the decaying form is essentially similar to the U = result. At the half-filling 
interband umklapp processes emerge and, according to the weak-coupling theory, open a 
charge gap, which results in an exponential decay of the pairing correlation. It is difficult to 
tell from our data whether P{r) decays exponentially. This is probably due to the smallness 
of the charge gap. In fact, Noack et al^ have obtained with DMRG an exponential decay 
for larger values of U, for which a larger charge gap is expected. 

When n is decreased down to 0.667 (Fig.^(b)), we again observe an absence of enhance- 
ment in P(r). This is again consistent with the weak-coupling theoryi: for this band filling, 
the number of electrons in the bonding band coincides with A^(= 30) at f/ = 0, i.e., the 
bonding band is half-filled. This will then give rise to intraband umklapp processes within 
the bonding band resulting in the 'C1S2' phase discussed in ref. ||, in which the spin gap is 
destroyed so that the pairing correlation will no longer decay slowlyll. 

In summary, we have seen that there are three possible causes that reduce the pairing 
correlation function in the Hubbard ladder: (i) the discreteness of the energy levels, (ii) 
reduction of Kp for large values of U/t, and (iii) effect of intra- and interband umklapp 
processes around specific band fillings. The first one is a finite-size effect, while the latter 
two are present in infinite systems as well. We can make a possible interpretation to the 
existing results in terms of these effects. For 60 electrons on 36 rungs with t± = 1.5t in ref. 



H| , for instance, the non- interacting energy levels have a significant offset ~ O.lbt between 
bonding and anti-bonding levels at the Fermi level, which may be the reason why the pairing 
correlation is not enhanced for U/t = 2. For a large U/t{= 8) in refs. p!o|JTT| , (ii) and/or (iii) 



6 



in the above may possibly be important in making the pairing correlation for t± = t not 
enhanced. The effect (iii) should be more serious for = t than for t± = 1.5t because the 
bonding band is closer to the half-filling in the former. On the other hand, the discreteness 
of the energy levels might exert some effect as well, since the non-interacting energy levels 
for a 32-rung ladder with 56 electrons {n = 0.875) in an open boundary condition have an 
offset of 0.15t at the Fermi level for t± = t while the offset is 0.03t for t± = 1.5t. 

Finally, let us comment on a possible relevance of the present result to the supercon- 
ductivity reported recently for a cuprate laddeiS, especially for the pressure dependence. 
The material is Sro.4Ca13.6Cu24O41.84, which contains layers consisting of two-leg ladders and 
those consisting of ID chains. Superconductivity is not observed in the ambient pressure, 
while it appears with Tq ~ lOK under the pressure of 3 GPa or 4.5 GPa, and finally dis- 
appears at a higher pressure of 6 GPa. This material is doped with holes with the total 
doping level of 6 = 0.25, where 6 is defined as the deviation of the density of electrons from 
the half-filling. It has been proposed that at ambient pressure the holes are mostly in the 
chains, while high pressures cause the carrier to transfer into the ladderJll. If this is the 
case, and if most of the holes are transferred to the ladders at 6 GPa, the experimental 
result is consistent with the present picture, since there is no enhancement of the pairing 
correlation for 6 = and 6 ~ 0.3 due to the umklapp processes as we have seen. Evidently, 
further investigation especially in the large-f/ regime is needed to justify this speculation. 

Numerical calculations were done on HITAC S3800/280 at the Computer Center of the 
University of Tokyo, and FACOM VPP 500/40 at the Supercomputer Center, Institute for 
Solid State Physics, University of Tokyo. For the former facility we thank Prof. Y. Kanada 
for a support in 'Project for Vectorized Supercomputing'. This work was also supported in 
part by Grant-in- Aid for Scientific Research from the Ministry of Education, Science, Sports 
and Culture of Japan. One of the authors (T.K.) acknowledges the Japan Society for the 
Promotion of Science for a fellowship. 
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FIGURES 

FIG. 1. The pairing correlation function, P{r), plotted against the real space distance r in a 
30-rung Hubbard ladder having 52 electrons for U = 1 with t± = 0.98 (□) and t± = 1.03 (O). The 
dashed line is the non-interacting result for the same system size, while the straight dashed line 
represents oc 1/r^. The solid line is a fit to the U = 1 result with t± = 0.98 (see text). The inset 

shows a schematic image of the discrete energy levels of bonding (0) and anti-bonding (tt) bands 
for f/ = 0. 

FIG. 2. A similar plot as in Fig.l for a 42-rung system having 76 electrons with t± = 0.99. 

FIG. 3. A similar plot as □ in Fig.l except U = 2 here. 

FIG. 4. The pairing correlation P(r) (□) against r for a 30-rung system fov U = 2 with (a) 
t± = 0.99 and 60 electrons (half-filled), and (b) tj^ = 1.01 and 40 electrons (half-filled bonding 
band). The dashed line represents the non-interacting result. 
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Fig.3 



10 



10 



-2 



10 




Fig.4(a) 
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Fig.4(b) 



